Lyapunov exponents and anomalous diffusion of a Lorentz 
gas with infinite horizon using approximate zeta functions 



Per Dahlqvist 
Mechanics Department 
^ ■ Royal Institute of Technology, S-100 44 Stockholm, Sweden 

On 

c 

a 
s 

(N 



O 



Abstract - We compute the Lyapunov exponent, generalized Lyapunov exponents 
and the diffusion constant for a Lorentz gas on a square lattice, thus having infinite 
horizon. Approximate zeta functions, written in terms of probabilities rather than 
] periodic orbits, are used in order to avoid the convergence problems of cycle expan- 

sions. The emphasis is on the relation between the analytic structure of the zeta 
function, where a branch cut plays an important role, and the asymptotic dynam- 
ics of the system. We find a diverging diffusion constant D(t) ~ logt and a phase 
transition for the generalized Lyapunov exponents. 
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^ ! 1 Introduction 

O 

Maybe the most well known measure of a chaotic system is the Lyapunov exponent. In the theory 
of chaotic dynamics one is of course interested in calculating this and similar quantities, either by 
finding analytical estimates or by devising effective calculation schemes but often one find oneself 
compelled to numerical simulation. This is unsatisfactory since it is not an easy task to extract 
information on the asymptotic behaviour from numerical data. 

Various averages of chaotic systems are obtainable via transfer operators and and their Fred- 
holm determinants or zeta functions. This is a beautiful formalism but the best results are 
obtained for a very restricted class of chaotic systems, namely those fulfilling Axiom-A. This 
is because Axiom-A guarantees nice analytical features of the zeta functions G|, ^, Q] which 
enables fast convergent cycle expansions to deduce its leading zeros Applications of cycle 
expansions to non Axiom-A systems are not very successful || 

In this paper we will study a system which is far from the textbook 1-d Axiom-A map, namely 
the two-dimensional Lorentz gas on a square lattice. This is a Hamiltonian system with two 
degrees of freedom, continous time and with an infinite symbolic dynamics. We will demonstrate 
that zeta functions may be of great use even here. The key point is that we will avoid writing 
the zeta function in terms of periodic orbits as that would lead us to divergence problems that 
we could not handle. The price we will pay is that our zeta functions are no longer exact, but 
they are approximate in a sense that won't effect the leading zero very much. The averages we 
will compute are directly related to the motion of this leading zero with respect to variations of a 
parameter. However, matters will be complicated if there are non-analyticities in the vicinity of 
the leading zero, like branch cuts. This is, we believe, a generic feature of non Axiom-A systems. 
Such singularities will cause problems for cycle expansions but do carry important information 
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about the asymptotic dynamics. In the Lorentz gas there will be a branch cut connected to 
the existence of the infinite horizon. It will not prevent the Lyapunov exponent from being well 
defined but will yield a diverging diffusion constant. It will also imply a phase transition for the 
generalized Lyapunov exponents. 

In section 2 we review the necessary theory. In section 3 we perform all calculations and we 
then end with some comments in section 4. 



2 Theory 

2.1 Lyapunov exponents and zeta functions 

The largest Lyapunov exponent is defined by 

A = lim -log|A(aso,t)| , (1) 

t—>oo t 

provided the limit exists and is independent of the initial point xq (except for a set of measure 
zero: namely the periodic points). A(xo,t) is the largest eigenvalue of the Jacobian along the 
trajectory starting at xo and evolving during time t. Using ergodicity this may be rewritten as a 
phase space average 

A = ~ f Li{dx)log\A(x ,t)\ = ~ <log|A(x ,i)| > , (2) 

where /i(dx) is the invariant density. 

We will now formulate the Lyapunov exponent in terms of evolution operators and zeta func- 
tions. Consider the following evolution operator acting on a phase space density according 
to 

= ( w(x,t)5(x-f(y)My)dy . (3) 



It is essential for the further development that the weight w(x, t) is multiplicative along the flow. 
We intend to evaluate the trace of this operator as that gives the expectation values of the weight 
(see appendix) 

< w >= lim trC* . (4) 

t— *oc 

The trace naturally turns out to be a sum over the periodic orbits 

trCl=J w(x, t)S(x - f(x))dx = E^E "g ^(1-mI) I ' (5) 

where n is the number of repetitions of primitive orbit p, having period T p , and M p is the Jacobian 
(transverse to the flow). w p is the weight associated with cycle p. 

Zeta functions arc introduced by observing that the trace may we written as 



For a Hamiltonian system with two degrees of freedom the zeta function reads || 

oo , —ikT \ m +l 



p 711 — 



where A p is the expanding eigenvalue of M p . 

Such infinite products over periodic orbits are often troublesome because they diverge beyond 
the first (nontrivial) zero. It is essential that the constant a in eq. (||) is sufficiently large so that 
the products converge. In our subsequent calculation it suffices if a is small and positive. 
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Cycle expansions of the zeta functions generally have better convergence properties since they 
converge up to the first singularity ||. The zeta function is entire for Axiom-A system which 
makes cycle expansions very successful for this special case. In this paper we will consider cases 
where the leading zero is also a singularity (branch point) so a cycle expansion will not converge 
even there and the zeta function is rather useless as it stands. We return to these problems in 
section 2.3. 

We must now find a weight w appropriate for computing the Lyapunov exponent. The 
quantity whose average we are going to study is Zo^|A(a;o, t)\ which is certainly not multiplicative. 
In one dimension we can study the average of the multiplicative weight w — |A(xq, t)\ T and obtain 
the Lyapunov exponent by differentiation 

1 d tr£} 

A=limi^lU. (8) 

t^oo t CIT 

The problem to find the appropriate weighted operator in more dimensions is nontrivial and 
the problem was recently solved |^|. But for our purposes it will suffice to insert the weight w 
above into the operator £} w just defined. The modification in || to make w exactly multiplicative 
is complicated but minor and we donjt expect it to affect the leading zero of the zeta function, 
and its vicinity, in which we are interested. 

The leading zero ko(r) = —i ■ h(r) is always on the negative imaginary axis, if r > 0. It 
provides the leading asymptotic behaviour of the trace provided that there is a gap until next 
zero or singularity. Then we have 

A=^(r). 0) 

Generalized Lyapunov exponents [jl0| A(t) are defined by considering the scaling behaviour of 
<|A(xo s t)r>: 

< |A(x ,t)| r >= lim tr4 = e A(r)rt , (10) 

t — »-oo 

so that 

A(r) = ^ , (11) 

T 

provided again that the leading zero is isolated. The ordinary Lyapunov exponent is recognized 
as the limit A = lim r ^o A(r). 

2.2 Diffusion coefficients and zeta functions 

We will consider the Lorentz gas obtained by unfolding the Sinai billiard. The coordinate in the 
unfolded system is called x. The corresponding vector in the billiard (or the unit cell) is x. They 
are related by translation x — x 6 T where T is the group of translations building up the Lorentz 
gas from the unit cell. 

The diffusive properties can be extracted from the average 

< e P< ft( - Xo) - Xo) > Xo . (12) 

The average is taken over one unit cell. Again we must perform the trick to introduce a multi- 
plicative weight and then by differentiation extract the average in which we are interested 

It was demonstrated in ref [fi"l[ that this average can be computed by considering the dynamics 
in the unit cell only. This is obtained by inserting the weight 

w(x,t) = e^K/V)-*) (13) 
into the evolution operator (^). The diffusion constant is now given by 

d = vt i, < (/ > o) - x ° )2 >-= & vt i, w tr ^ • (l4) 

i=i i=i ^ i 



where the sum is taken over spatial components of the 2v dimensional phase space. 
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2.3 Approximate zeta functions 

In some recent publications we have investigated a way of approximating zeta functions for 



intermittent systems |7j, |12|, 13 . We call it the BER- approximation after the authors of rcf 
14 . In an intermittent system laminar intervals are interrupted by chaotic outbursts. Let A$ 
be the time elapsed between two successive entries into the laminar phase. The index i labels 
the i'th interval. Provided the chaotic phase is chaotic enough, the lengths of the intervals A; 
are presumed uncorrelated, and A may be considered as a stochastic variable with probability 
distribution p(A). The zeta functions (unit weight w = 1) may then be expressed in terms of the 
Fourier transform of p( A) 

/>oo 

Z(k) sa Z(k) =1- e~ ikA p{A)dA . (15) 



Due to the normalization of p(A) we see that leading zero kg = is by construction exact because 
of probability conservation. 

In order to compute the probability distribution we introduce a surface of section (SOS). This 
should, according to the BER prescription be put on the border between the laminar and chaotic 
phase. We call the phase space of the SOS f2 and its coordinates x s . The flight time to the next 
intersection is then a function of x s : A s (x s ). The probability distribution then reads 

p(A) = f 6(A-A s (x s ))fM(dx s ) , (16) 
Jn 

where fj,(dx s ) is the invariant density, which is uniform n(dx s ) = dx s / J^dxg for Hamiltonian 
ergodic system, assuming of course that the SOS coordinates are canonically conjugate. 

It is straightforward to include weights in this formalism, like e.g. w = \A(xo)\ T . The local 
expansion factor A s (x s ) is also a function of x s . The zeta function is then related to a generalized 
distribution 

p T (A)= [ \A s (x s )\ T S(A-A s (x s ))dx s . (17) 
Jn 

The zeta function Z T (k) is obtained by inserting p T (A) into jig). 



3 Application to the Lorentz gas 

We now begin our study of the Lorentz gas on a square lattice. The lattice spacing is unity, the 
disk circular with radius R < 1/2, and the point particle bouncing around in this array has unit 
velocity. It is the unit cell of this system, the Sinai billiard ]l5[ whose dynamics we will study and 
this is indeed an intermittent system; there exist periodic orbits with arbitrary small Lyapunov 
exponents logA p /T p . The disk will define our SOS. We use the two angles <fi and a defined in 
fig la as coordinates. The normalized measure is then dx s = d<p d(sina)/47r. Consider now a 
segment of the trajectory between two disk collisions. This segment can be labeled according to 
the disk that would be hit in the unfolded system, the label q = (n x ,n y ) is the associated lattice 
vector ]l2] ]. It is easy to realize that only disks associated with coprime lattice vector may occur. 

The purpose is now to apply the BER approximation to this system and compute Lyapunov 
exponents and diffusion constants, cf refs |l^, 



3.1 Calculation of p r (A) 

In this section we will derive the following approximate expression for p T (A) 

Pt(A)~\ 4 r( 2_, )2 2T/2 _ 3fi _ T/2 _ 2 ^> 

3 r(2-r) A 3 -3-/ 2 A>l/^Yt . 
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Figure 1: a) The Sinai billiard with definitions of the variables 4> and a. b) The unfolded 
system with free directions (corridors) indicated, c) The region of integration. 
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This expression was used already in ref p3[ but as it plays a central role in this paper we present 
its derivation in some detail. Complementary details may be found in refs jij], |l3| . 

We start from expression (|l7| ) for p T {A). First we partition the SOS into subsets f2 q where 
fi q is the part of f2 for which the trajectory hits disk q. Then we smear the distribution, that is, 
we replace the delta function with some extended distribution S a . The exact form of this function 
is irrelevant, the only thing we assume is that the width is big: a ^ 1. We will be interested in 
the behaviour of the leading zero of Z T (k) for small r. This zero is then close to the origin and 
it is evident that smearing of p T {A) will only have minor effect on this zero. 

We now have 

p T (A)=Y / [ \A s (x s )\ T 6 a (A-A s (x s ))dx s . (19) 
q Jn« 

The large width a allow us to move the smeared delta function to the left of the integral sign 
because the variation of A s (a; s ) over il q is of the order ~ R and we have a ^> 1 > R. 

p T (A)=J2^-l) [ \As(x s )\ T dx s =^<5 CT (A-g)a q (r) . (20) 

q q 

We have chosen the length of the lattice vector |q| = q as a mean value of A s (x s ). The local 
expansion factor is |A s (cc s )| = 2A s (x s ) / Rcos(a) . 

It is easily shown that the phase space area taken up by disk q is given by the inequality 

| ^sin(</)-0 q -a)+sin(a;) |< 1 , (21) 
R 

where # q is the polar angle of the lattice vector q. Generally parts of this region are eclipsed by 
disks closer to the origin. So in order to find f2 q one has to subtract these. We will focus on the 
limit of small R which gives the more easily handled inequality 

||(0-0 q -a)+sin(a)| < 1 . (22) 

Let us begin with the limit of small A. In ref |l2| it is shown that if disk q lies within a 
certain radius: q < 1/2R, there are no eclipsing disks in front of it, and expression ( p2[ ) may be 
used directly, and the integral is easily calculated: 

In order to find an approximate expression for p T (A) we must know the density of coprime lattice 
points, i.e. the average number d c (r)dr of such points having a distance between r and r + dr 
from the origin. In ref jl3j this was, to leading order, found to be d c (r) ~ This yields 

p T (A)=^^(A- (? ) aq (r) = g^^i? 1 -W . (24) 

q ^ > 

Obviously the disk radius has to be small for this to apply. 

Next we consider the opposite limit A ^> 1/2R. For each (coprime) disk q fulfilling q < 1/2R 
there are two (or one, depending on symmetry) transparent corridor in the direction q |]l2| , |l6)| , sec 
fig lb. Far beyond this critical radius the accessible disks will be those adjacent to the corridors. 
(They still have to coprime so they will lie on one side of the corridor only, see fig lb). We will 
discover that this will lead to a power law decay of p(A). We will be interested in the particular 
power (as a function of r) and not the prefactor. For that reason we perform our calculation 
in the corridor having direction vector (1, 0) as all corridors provide the same power. In this 
corridor the accessible disks are the ones being labeled (n, 1). Disk (n, 1) is shadowed by (1,0) 
and (n — 1, 1). We need to evaluate the integral (cf eq (|23|)) 

cos T+1 (a)da d<f) . (25) 
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' Numerical' 
Theory 



Figure 2: Lyapunov exponent versus disk radius according to numerical simulation and 
eq (32). 



It is more convenient to consider the sum 



j n = J2ii = J 



cos r+1 (a) da d<fr 



(26) 



because this integral has support from the triangular region in fig I c (it is a triangle only in the 
limit n — > oo of course). 

From eq (^2|) one can deduce that the base length (in the <j> direction) of this triangle scales 
as ~ 1/n and the height ( in the a direction) as 1/y/Jn). A short calculation now yields that 
J„ ~ l/n 2 ~ r / 2 . Differentiation gives j„ ~ 1/n 3-7 "/ 2 . The fact that Z q ~ n together with eq. f ]23| ) 
implies that the a q 's decay as ~ l/n 3 ~ 3r / 2 . The density of accessible disks is uniform in A (since 
they lie along the corridors) so our final result is 



P (A) 



A3-3T/2 



A » JR. 



(27) 



This power law does not depend on the small R limit. 

The prefactor of eq (|27|) is determined by demanding that _p T (A) is continous at A = 1/2R. 
In order to get correct normalization for r = we multiply the entire p T (A) by 13/12. This 
approximation means a rather crude neglect of the of the transitional behaviour at A « 1/2R. 
In the following we must be very cautious when using it since some results may be more sensible 
to our approximation than other. 



3.2 Lyapunov exponents 

The approximate zeta function we are going to work with is obtained from (|1 



z T (t)« 1-/; 



-ikA 



MA) 



3 r(2— r) "2 T + 



r (« -1 - t 7(t + 1,z) + z 2 - 3r / 2 r(2 + 3t/2, z)) , 



(28) 
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where z = The functions 7(0, z) and T(a, z) are incomplete gamma functions [fL7| . Expanding 
this to first order in r gives 

Z= {z + 4(logz + 7 -^)...} + 

{(-^+log(2 J R 2 )) + (f +log(2 J R 2 ))z...}r ... 

The derivation of this expansion is rather lengthy but the only step that is slightly tricky is 
the expansion of the incomplete gamma function T(a, z) near an integral power a. We have for 
convenience thrown away all indices of the zeta function and use strict equality sign but we must 
not forget that we work with an approximation of an approximation of the exact zeta function. 

We must now proceed with some care since the leading zero is not isolated, a branch cut 
along the negative real z-axis reaches all the way up to it (we have choosen the principal branch 
of the logarithm) . 

We are interested in the following derivative of the trace, cf eq (|8|) 

d ■ - 1 _ d 1 f *t' d ,„„ ry , 1 1 f nZt > d d 



trC 1 = ^~^ e zt —log Zdz | T=Q = — / e zt ^^logZ | T=0 dz , (30) 
dr dr Zm J dz Zm J dr dz 

where we have differentiated inside the integral sign. We have now formulated eq (^) in the 
rescaled and rotated z — plane using rescaled time t' = 2Rt. The function to be Fourier trans- 
formed is 

dr dz 12 z + ^-(logz + 7 - -^-) 2 . . . 12 z z 3z 

(31) 

where we have kept only those terms yielding the leading term and the first correction. Collecting 
it all together yields 

A = lim2i?(^ - log(2i? 2 ))(l + ...) = 2R(^ - log(2i? 2 )) . (32) 

The limiting value is due to the behaviour of the zero but the power law correction are due to the 
fact that it sits on a branch point. The particular size of the first correction is not very accurate, 
as we will realize after reading section 3.4 and 3.5, but the important thing to bear in mind is that 
there exist slowly decaying corrections indicating slow convergence of the Lyapunov exponent in 
numerical computations. 

In fig 2 we compare numerical results on the Lyapunov exponent with our expression A = 
2R(^ - log(2i? 2 )). The numerical values, from (jj), 

are calculated for rather large disk radii, 
where we should not expect much agreement a priori. Nevertheless our estimate is only 5% wrong 
when R = 0.1. The reason why the the numerical values exceed our estimates for large R is easily 
understood. This is because the disk faces (in the unfolded system) come closer to each other so 
that taking the length of the relevant lattice vector (as we did) overestimates the time of flight 
between them. 

The Lyapunov exponent of the corresponding Poincare map with the disk defining the SOS 
is related to the Lyapunov exponent of the flow according to (l9) \ map — A < A s A/2i? w 
(s ~ l°g(2-R 2 )) where < A s > is computed by means of expression (|l8|). We see that A — > 
— 21og(i?) + 7/12 — log2 when R — > which agrees with the conjectured limit ]2(], EH]]: A — > 
— 21og(i?) + C + O(R). Indeed we have analytically found an estimate of the constant C » 
7/12 — log 2 « —0.110. This is very close to the numerical value found by p0| |, as far as we can 
extract it from fig 2 in ref p0[). 

In fig 3 we plot the generalized Lyapunov exponents for different disk radii. Note that when 
r > the leading zero is indeed isolated and we do not have to worry about the cut. The position 
of the zero is computed numerically. A(l) is the topological entropy which tend to a finite limit 



when R — * 12 , whereas A(0) — > as R — > 0. When r < the branch cut itself will provide 



the leading behaviour of the trace - a power law 12, 13], and the generalized Lyapunov exponent 



will be zero. This means that A(r) cannot be analytic at r = 0. This is referred to as a phase 
transition Jl(], . 
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Figure 3: Generalized Lyapunov exponents A(r) versus r for two disk radii. 



3.3 A first calculation of the diffusion constant 

In order to study diffusion we calculate the generalized probability distribution p w (A) using 
appropriate weight (^). However, we only keep the spatial components of (3 yielding the two 
dimensional vector j3. As we study smeared p^(A) it suffices to approximate the spatial part 

of (/*(a;) — x) with the lattice vector q. So, we must now compute the generalized probability 
distribution p^(A) a la section 3.1 but using the weight exp(/3 ■ q) = ex^{[3qcos{4>fj — 4> q ) where 
we have written q = g(cos(0 g ), sin(</> g )) and /3 = /3(cos(^), sin(<^)). We can use the result 
of section 3.1 to some extent, since we realize, after inspecting eq (|T^), that we can make the 
following factorization 



P/3 (A) =po(A) 



I 



A/3eos( 



(33) 



fd(f> q 

To know the support of this integral we need to know about the angular distribution of access- 
able disks for a particular value of q = A. To this end we introduce the following simplifying 
assumptions 

1. the coprime lattice points are distributed isotropically. 

2. if q > 1/2R we assume only four corridors with 4> q equal to one of 0, 7r/2, tt or 37r/2. 

Assumption 2 means a rather brute neglect of the transitional behaviour at A w 1/2R. We 
neglect the important fact the number of corridors grows when R — * 0. We comment on this in 
section 3.5. The important thing now is that we preserve the symmetry in order to prevent net 
drift. 

We first consider the case q < 1/2R. Using assumption 1 above gives 



t R ±. f V Acos 
3 2vri 



= -RIo(PA) 



-i?(l 



/3 2 A 2 



•) 



(34) 
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where Iq is a modified Bessel function In order to calculate the zeta function we need the Fourier 
transform of this 

l/2fl r, on 

Pf3 (A)e~ lkA dA = -(1 - z/2 + z 2 /6 . . .) + ^-(1/3 - z/4 + z 2 /10 • • ■) , ■ ■ ■ (35) 

where we have used rescaled variables z = ik/2R and /3' = (3/2R. 

Next we consider the limit q > 1/2R. We find, using the second assumption above 

P /3 (A) = 3^i ( ^ A + e ~" LA + ^ + e ~ ' yA) A > « ■ (36) 
Fourier transforming this gives 

P0 (A)e~ ikA dA = i(£ 3 (z - /£) + E 3 (z + X ) + E 3 (z - /? ) + E 3 (z + (3' y j) . (37) 

1/2R O 

The resulting zeta function can now be expanded 



Z (3 = 1 - O^AK^dA = i( 7 - ^)/3' 2 + z(l - m + z 2 I( 7 - 

i [(z - /3^) 2 log(z - + (z + p x f \og(z + 0> x )+ 

(z - /?;) 2 io g (z - p v ) + (z + /?;) 2 io g (z + /?;)] 



(38) 



According to sec 2.2, we are interested in the following quantity 

(|j + |*)ir£ 1^=0= + |^ / A log Z dz l^o 



(39) 



We now proceed in the same way as in sec 3.2. We want to determine the asymptotic behaviour 
of the fourier transform of 

( ^ + ^ ) ^ l0gZ| ^i^- 1 ° gZ - 7) ' (40) 
To this end we now need the following integral 



1 f logy.- 
2iri 



dz = t'(l - 7 - logt') . (41) 



In deriving this it is convenient to use contour integration and let the contour encircle the negative 
real z axis. We thus find the following diverging diffusion constant (D = limt_ (00 D(t)) 

m = k { Wi + M )trC b '=°= + 3,09(2R!)) ' <42) 

This logarithmic divergence of the diffusion constant agrees with ref |l6| but the prefactor is 
not correct. The computation will be refined in the section 3.5. The important thing to learn 
from this section is that the Laplacian exposes the logarithm in the series expansion of the zeta 
function. So we now know where to focus our attention, namely at the tail of pg(A). 

3.4 A closer look at the tail 

At the end of section 3.1 we aimed at finding a good approximation for p r (A) for the whole range 
l/2i? < A < oo. In this section we will make a more careful investigation of the A — ► oo tail of 
p T (A). We restrict ourselves to the case r = and will compute the limit limA^oo A 3 po(A) for 
small disk radii. 
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Let us look at the corridor with direction vector q = (n x , n y ) where n x and n y are coprime. 
Suppose for a moment that q lies in the first octant so that n x and n y positive and n x > n y . The 
accessible disks in this corridor are the ones labeled q' + nq and q" + nq where q' and q" are the 
predecessor and successor in the Farey sequence of order n x |T^j , A calculation analogous to the 
one at the end of section 3.1 now gives the following expression for a q (0) (as defined in eq (|23|)) 

1 2qR + ^ - 2 



--« = t ' /S +0(1/n4) (4:J) 



and the same holds for the sequence q" + nq. The contribution to po(A) is 

V+n^(A - W + nq\) = ^3 ^( 2 ^ + " 2 ) ( 44 ) 

n 

To compute the tail of Po(A) we need to sum over all coprime q such that q < 1/2R. We restrict 
the summation to the first quadrant q G S 

S = {q = {n x ,n y )\n y > 0; n x > 0; (n x ,n y ) = 0; ^Jn 2 x + n 2 v < 1/2R} (45) 

and multiply the result by four, to account for all four quadrants, and then by two, to account 
for both q' and q" defined above. The result is 



MA) - ^ J>Jl + ^ - 2) + 0(1/A 4 ) (46) 



The nontrivial sum above is treated by Bleher [[16| and we insert the small R limit of the sum 
and obtain 

MA) ^ A52^ +0(1/A4) (4?) 
Suppose that we expand the unweighted zeta function Zq into the more general series 

oo oc 

Z = Y,^z l + Y,^ z ^ogz (48) 

i=l i=2 

For the approximation ofp(A) worked out in section 3.1 we have b{ — for i > 3. The correct value 
of &2 may be computed from eq (|47]) (using e.g. standard expansions of exponential integrals) and 
is found to be 62 = l/V 2 which differs considerably from the result in section 3.1. The reason for 
this error is the absence of a smooth transition at A = 1/2R in eq ([T8|). The crossover behaviour 
is slow as indicated from the (9(1/A 4 ) term above. Generally we expect the asymptotics of the 
tail to look like p(A) ~ J2 n >3 °nl A n so that there now may exist nonzero bi of any order. 

3.5 A second calculation of the diffusion constant 

One could suspect that the coarse assumptions, especially assumption 2 in section 3.3, induce an 
error in the prefactor of super diffusion. We will now demonstrate that the prefactor is robust 
against refinements of the this assumption. A better estimate of pp(A) is obtained by replacing 
exp{(3' x A) = exp(/3'A cos(<pp)) in eq (J3q) with the averaged quantity 



1 



2c/A 



c/A 
-c/A 



The integration range decreases like 1/A since the width of the corridor is constant. For suffi- 
ciently large A we can take j3' cos{<f>p — <f>) r* P' x + (3 y (f> and find that the integral above is 



e^ A (l + 0((3' y 2 ) 



(50) 
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The correction is 0((3' 2 ) and has no effect on D(t), as we anticipated. 

We also reduced the number of corridors to four, irrespective of the disk radius. As long as we 
are not interested in the angular properties this is easily justified by studying how the laplacian 
acts on (|38|). 

We are led to the conclusion that the prefactor of super diffusion depends solely on the 
coefficient hi of the unweighted zeta function ( |i"8| ) according to D(t) = ^\ogt. Inserting the 
correct value of 62 from the previous section we arrive at 

m = 2]h logt (51) 

which is indeed the exact result 

Again there are slowly decaying corrections to D{t) and it is not surprising that the numerical 
detection of this diffusion law has eluded serious attempts [^U . 



4 Concluding remarks 

The main advantage of the approximation outlined in this paper is its apparent simplicity. It is 
far more easy to apply than a proper cycle expansion would be, taking all the (infinite number of) 
pruning rules into account, and it gives a very good description of the leading zero. This zero and 
its vicinity, was our only interest in this paper and we could work out a very coarse description 
ofp(A). 

In ref. we demonstrated that refinement of the expression for p(A) down to a scale ~ R 
also accounted for a good description of the trace down to this scale. This means that some 
nonleading zeros may also be extracted within the BER approximation. However, in order to 
compute higher order spectra, like semiclassical spectra, one needs to take the correlation between 
laminar intervals into account. The way to do this would perhaps be to combine the knowledge 
of the asymptotics on the periodic orbits, as obtained from the present approximation, with a set 
of explicitly computed periodic orbits. 

The major drawback of the BER approximation is the lack of error term; so far we have 
no bounds on the errors in our computations. It is of course highly desirable to work out such 
bounds which would put the present considerations on mathematically firm ground. 
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Appendix 

In this appendix we will verify equation (Q) < w >= limj—Kx, trC l w . 

The invariant density may be expressed in terms of periodic orbits as fi(dx) — lim^oo /j, t (dx) 
where 



S(t — nTp) f f ,^dx' 
- \det(l-M?)\ 

p^p.p.o. n—1 t 



= E E iJn-J^i L/^-^v (52) 



where v is the speed at x. 

The invariant density is a eigenfunction of the evolution operator (w = 1) corresponding to 
the leading eigenvalue (equal to one). This requirement is satisfied by /it(dx) for any t. However, 
/it(dx) is only a smooth function in the limit t — > oo (we refrain from showing this). It follows from 
the boundedness that the density is normalized J fit(dx) — * 1 as t — ► oo [?, |||. The expectation 
value can now be expressed in terms of periodic orbits as 



< w >— j w{x,t)^{dx) = (53) 
6[t-nT p ) 



Hm y r p r^" ,;v ;\, (54) 

Gp.p.o. n—1 y 



Comparing this expression with eq (|5|) immediately verifies eq 



